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ABSTRACT 


In this paper, the voltage stability of the power system is studied during fault 
conditions. Enhancement of the system’s stability will be achieved by 
utilizing Flexible AC transmission systems (FACTS) controllers at the best 
place in the system. The optimum placement of FACTS controllers occurred 
on the most affected bus by the fault in the system (weakest bus). Two 
approaches have been used in this study to effectively obtain the best 
location of the FACTS controllers in the system. The first method is based on 
computing the deviation which occurs in the active power and reactive power 
due to the fault at each load bus at a time. Whereas the second method is 
performed through exploiting a MATLAB fuzzy set technique utilizing two 
indices: Line Flow Index (LFI) and Voltage Profile Index (VPI) during fault 
and steady-state conditions. The results show that both of these indices 


STATCOM resulted in the same bus as the best location. Remedial actions in the attempt 
SVC at improving in the stability of the power system are discussed taking 
VPI the advantage of using FACTS compensation SVC and STATCOM at the 
most vulnerable system buses. In this work, MATLAB program with an 

IEEE 24 bus system is examined. 
This is an open access article under the CC BY-SA license. 
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1. INTRODUCTION 

In current years, due to the fact that increasing of power demand and limiting the expansion of 
the power generating units and transmission lines because of environmental and economic reasons, several 
transmission lines are reaching to their limits and the power systems enter to instability problems. In 
principle, the electrical system must deliver trustworthy power to the consumer, therefore the reliability of 
the system must be improved by transferring the power to places of use. Numerous developed methods have 
been launched in terms of a support system. As a result, system stability becomes a very challenging problem 
that must be considered [1-3]. The voltage stability analysis is crucial issues because unstable voltage leads 
to the complete or partial interruption of the power system. Also, the Reactive power plays a significant role 
in the constancy of the power system [4, 5]. The failure of the power system to encounter the demands 
reactive power in the high load systems, it is one of the main reasons for the instability of voltages in the 
system. The only way to solve these problems is to decrease reactive power load or increase extra reactive 
power before reaching to collapse point [6, 7]. Recently with the evolution in control theory, in addition to 
the electronic revolution witnessed by the world, which have allowed the design and implementation modern 
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controllers known as Flexible A.C Transmission System (FACTS). The FACTS controllers are primarily 
employed for solving numerous power system control problems. However, modern studies show that the 
possibilities of employing FACTS controllers to develop power system stability in supplement to their core 
function of power flow control. Also, FACTS devices could be used to improve voltage profile, power 
transmission capacity, and enhancing power system stability which is a suitable control strategy to enhance 
different system performances such oscillations, losses of the system, stability, reliability, security, and 
power factor in distribution systems [1, 8-10]. Although the FACTS controllers are expensive, they provide a 
and fast and smooth response to provide stability for the power system during the transient and steady-state 
and allow utilities to run their transmission schemes even nearer to their physical limits, compared with 
conventional methods [11-13]. The location of FACTS controllers on the power system are very important 
factor during stability improvement of the system. There are numerous techniques for locating FACTS 
devices such as contingency sensitivity index (CSI), Genetic Algorithm (GA), sensitivity index, line outage 
distribution factor (LODF), and Particle Swarm Optimization (PSO) [14-16]. In order to improve the voltage 
profile in the buses- reactive and active power under diverse loading circumstances-, FACTS devices with a 
fuzzy system can be used as static controller which leads to advance the stability in power electrical 
systems [17]. To achieve the stability, security requirement, controllability and maximum system loadability 
in a transmission line of the system are achieved by an optimal logic fuzzy PID, and genetic algorithm (GA) 
controller is designed for Thyristor based FACTS controller [18]. In [19], Genetic Algorithm (GA) is utilized 
to find the optimal location of FACTS device based on MATLAB program for satisfying the stability of the 
system. In [20, 21], Gravitational search algorithm (GSA) optimization methods is presented for selecting the 
optimal location. 

The aim of this paper is to select the best location of the FACTS controllers on the power system to 
expand the performance and stability of the system during the steady state and transient state. There are two 
methods will be used to get the optimum location of the FACTS: the first method is based on the computing 
the deviation occurs on the active power and reactive power at load buses due to the fault at each load bus at 
the same time, where the bus has the highest deviation will be considered the weakest bus while the bus has 
least deviation is considered the most robust bus. The bus that is most affected during the fault and has the 
highest deviation is considered the weakest bus and the optimum placement for FACTS controller in the 
system. The first method is called the traditional method, which is occur only during transient state condition 
(fault). However, the second method is intelligent method, achieved by using a MATLAB fuzzy set 
technique which are used two indices: Voltage Profile Index (VPI and Line Flow Index (LFI). It depends on 
the value of CI index, where the bus which has the highest CI value is ranking the first and considered the 
weakest bus, also it is considered the best location of the FACTS devices in the system. This method can be 
used during steady-state and transient state (fault) as will be described in the next section. 


2. FACTS ALLOCATION METHODOLGIES 
2.1. Traditional optimum FACTS allocation methodology 

This method is called a traditional method. It is based on computing the deviation occurs in the 
parameters of each load bus (active power and reactive power) due to the fault at each load bus one at a time. 
The deviation of these parameters can be computed by (1) and (2) below, where the load bus that has the 
highest deviation will be considered the weakest bus and the perfect place for FACTS devices. This method 
can be used during the fault condition only. 


AP = Prout — Poase (1) 
AQ T Q fault E Qbase (2) 


where: 

AP, AQ are active power and reactive power deviation respectively. 

Prautt, Qfaut are active power in MW and reactive power in MVar during the fault. 

Phase, Qbase are active power in MW and reactive power in MVar pre fault during steady-state condition 


2.2. Fuzzy optimum FACTS allocation methodology 

A fuzzy logic controller is influential tool for regulating complex processes and able to clarify the 
complicated issues without the necessity of a detailed statistical model of the system. The durability of the 
varying network specifications is another merit of the fuzzy logic controller [22, 23]. In this method, they are 
adopted two indices LFI and VPI, that correspond to the margin of voltage stability for a selected operational 
condition. The initial index is the “Line Flow Index” (LFI) that identified by (3) below, the derivation at: 
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LF] = 2P:l1+cos(@-®)] (3) 


Vs? YLcos(®) 


When (LFI) is near to unity, that means the maximum loadability point is reached, after this point, 
the system becomes critical [24-26]. By utilizing fuzzy setting within the MATLAB model program, LFI is 
classified into five curves (memberships) to characterize the input of the fuzzy laws as shown in Figure 1. 
The output rules of fuzzy are representing the impact of line flow either due to a high load during steady-state 
or due to fault during the transient state, and it is known as "Severity Index of Line Flow" (SILF). SILF is 
also classified into five curves (memberships) as shown in Figure 2. The voltage in p.u. at load buses of the 
system is the second index, and it is described “Voltage Profile Index” (VPI) as shown in Figure 3. The VPI 
supplied to the rules of fuzzy as four input curves (memberships) to give the fuzzy output which call 
"Severity Index of Voltage Profiles" (SIVP). It is classified into five curves (memberships) as shown in 
Figure 4. Can be determined the most vulnerable system buses by Criticality Index (CI), which computes by 
summation the two severity indices SILF and SIVP. According to CI value will be able to determine the 
weakest bus then ranking all load buses. Where the bus which has the highest CI value is ranked the first and 
considered it the weakest bus. In the Figures 1- 4, the symbols that mean: Small index (S), Very Small index 
(VS), Medium index (M), Very High index (VH), High index (H), Low Severity (LS),Very Low Severity 
(VLS), Above Severity (AS), Below Severity (BS), More Severity (MS), Normal Voltage (NV), Medium 
voltage (MED), Low Voltage (LV), Over Voltage (OV), PRE-Severity (PRES), More Severity (MOS) and 
Medium Severity (MEDS). 
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Figure 1. LFI curves Figure 2. SILF curves 
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Figure 3. VPI curves Figure 4. SIVP curves 


The proposed method is explained in the next steps: 
— Selecting load bus one by one (one at a time), reactive power demand is increased until LFI reaches to 
the extreme value. 
— Calculating LFI and VPI. 
— LFI and VPI are fed in fuzzy set notation to give SILF and SIVP by using fuzzy rules, then calculate Cl 
for each load bus by (4) below. 


CI=¥ SILF + ¥ SIVP (4) 


— Repeating the above steps for each load bus. 
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— Ranking all the buses descending based on the value of the CI for each load bus, where the most 
significant value of CI is the bus being the weakest. 

— During the fault conditions, the steps from step2 to step5 will be repeated. By assuming the fault 
influence on load buses is similar to the severity which occurs on load buses during maximum 
loadability. 


3. RESULTS AND DISCUSSION 

In this paper, IEEE 24 bus system will be used. The buses of the system are connected by 
transmission lines. Bus13 is a slack bus, buses (1, 2, 7, 14, 15, 16, 18, 21, 22 and 23) are PV buses and buses 
(3, 4, 5, 6, 8, 9, 10, 19 and 20) are load buses. The system has 11 generations, 9 loads; each one has active 
MW and reactive power MVar and 38 transmission lines as explained in Figure 5. 
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Figure 5. Single line diagram of IEEE 24-bus system 


3.1. Steady-state results 

MATLAB simulation results during steady-state will be demonstrated in Table 1 below. These 
results represent Voltage (Vin p.u), Active power and Reactive power values at each bus in the system during 
steady state. To determine the optimum placement of the FACTS devices, the intelligent method will be 
used, achieved by using a MATLAB fuzzy set technique where, LFI and VPI will be calculated at each load 
bus, one at a time, during the severest loading condition (high load point). Then, the two inputs LFI and VPI, 
for each load bus, are fed into fuzzy rules to give SILF and SIVP respectively. By summing SILF and SIVP 
will be getting to CI for each load bus during high load point to control the weakest bus and rank the load 
buses according to CI value from high to low value as shown in Table 2. 

Table 2 is demonstrated the CI value of each load bus, then ranking the load buses based on its 
value. Therefore Bus 9 is ranked at first and it is considered the weakest bus, because it has the biggest CI 
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value. Then, it is selecting as the suitable placement for SVC and STATCOM device to enhance the stability 
of the system 


Table 1. Steady-state load flow Table 2. CI and ranking during steady-state 
Busno. P[MW] Q[Mvar] Voltage [p.u] Bus no. CI=)'SIVP + )'SILF Rank 
1 -89.27 -766 1.035 9 1034.2 1 
2 -99,4 -195.1 1.035 3 1032.09 2 
3 204.3 -241.9 0.9281 6 1022.79 3 
4 36.25 88.15 0.9439 10 1010.18 4 
5 -12.9 -164.4 0.94 5 1003.25 5 
6 -78.45 -44.82 0.8644 4 986.44 6 
7 63.55 342.4 1.025 8 939.81 7 
8 -404.4 -53.56 0.9448 19 906.14 8 
9 -298.1 344.3 0.8973 20 898.13 9 
10 -291.2 293 0.8513 
11 -289.4 76.26 1.01 
12 -301.2 66.3 1.03 
13 197.2 -14.32 1.02 
14 378.9 -24.24 0.9881 
15 -225.8 -124.1 1.014 
16 188.3 -43.93 1.017 
17 -338.6 -81.08 1.035 
18 -62.93 -104.4 1.05 
19 -182.3 -19.25 1.024 
20 -137.5 -18.64 1.039 
21 -399.9 -124.8 1.05 
22 -299.9 21.86 1.05 
23 -364.6 -110.4 1.05 

24 -204.6 8.668 1 


3.2. Results during fault (Transient State) 

Applying three phase faults at system load buses, one at a time, during 0.3 sec without inserting any 
compensation components. During these conditions, will be flow the two methods: intelligent MATLAB 
fuzzy method and traditional method to determine the optimum location of FACTS devices. Now by fuzzy 
methodology the LFI and VPI will be also calculating, and they are also fed into fuzzy rules to give fuzzy 
outputs indices SILF and SIVP respectively. Summation of this indices will be resulted CI for each load bus, 
the load buses will be ranked depending on the CI value for each load bus, and the bus which has the highest 
value considered the best location of the FACTS devices and the weakest bus in the system as demonstrated 
in Table 3. 

From the obtained results, as shown in Table 3, Bus 9 is ranked at first because it has the highest CI 
value, therefore it considered the weakest bus, whereas Bus 20 is ranked at last because it has the lowest CI 
value, therefore, it is considered strongest one. Therefore, Bus 9 is the optimum location of SVC and 
STATCOM controllers. Based on traditional method during fault occurred at load buses respectively, then 
computing the deviation value occurs in the parameters of each load bus (active and reactive power) by using 
(1) and (2) of AP, AQ at Bus 9 has highest deviation, while Bus 20 has a value of the least deviation as shown 
in Figure 6 , Figure 7 and Table 4 , Table 5 below. That means Bus 9 is the most vulnerable bus and it is the 
best placement for SVC and STATCOM controllers. 


Table 3. CI and rank during fault condition 


Bus no. CI=)'SIVP + )'SILF Rank 
9 919.86 1 
3 919.65 2 
10 918.29 3 
19 918.12 4 
6 918.09 5 
4 917.96 6 
5 917.47 7 
8 917.24 8 
20 906.14 9 
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From the results, as shown in Figure 6 and Table 4, noted the highest active power Deviation (AP) 
occur in Bus 9, while Bus 20 has a value of the least active power Deviation AP. That means Bus 9 is the 
most vulnerable bus and it is the best placement for SVC and STATCOM devices. 
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Figure 6. AP in MW due to fault at load buses respectively 
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Table 4. AP in MW due to fault at load buses respectively 
Bus no. Active power deviation (AP = Pfault — Pbase) [MW] during fault at: 


Bus3 Bus4 Bus5 Bus6 Bus8 Bus9 Bus10 


3 +2.2 +1.7 +1.5 +1.8 +0.6 +2.3 +1.8 
4 -0.58 -0.47 -0.37 -0.49 +041 -0.62 -0.46 
5 +0.87 +0.7 +0.53 +0.72 +0.48 +0.93 +0.69 
6 +0.06 +0.06 +0.03 +0.04 +0.21 +0.07 +0.04 
8 -0.4 -0.3 -0.3 -0.3 +0.9 +0.4 +0.4 
9 -4.4 -3.3 -2.8 -3.5 -1.2 -4.6 -3.6 
10 -3.9 -2.9 -2.5 -3.1 -1.1 -4.1 +3.2 
19 0 0 +0.1 0 +0.1 0 +0.1 
20 0 0 0 0 0 0 0 


Bus19 
+2.7 
-0.64 
+1.05 
+0.09 
-1.3 
-5.3 
-4.7 
-0.2 
-0.1 





Bus20 
-7.2 
+2.33 
-3.47 
-0.26 
+0.8 
+15.6 
+13.7 
-0.1 
+0.1 


From the results, as shown in Figure 7 and Table 5, noted the highest reactive power Deviation AQ 
occur in Bus 9, while Bus 20 has a value of the least reactive power Deviation AQ. That means Bus 9 is the 
most vulnerable bus and it is the best placement for SVC and STATCOM devices. Based on all result, Bus 9 
is the best location for SVC and STATCOM devices during steady state and transient state which achieved 


by intelligent MATLAB fuzzy method and traditional method. 
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Figure 7. AQ in MVar due to fault at load buses respectively 
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Table 5. AQ in MVar due to fault at load buses respectively Steady-state load flow 


Bus no. 


CONN BB W 


19 
20 


Reactive power deviation (AQ = Qfault — Qbase) [MVar] during fault at: 


Bus3 
-3.1 
+1.19 
-1.8 
+0.03 
-0.06 
+3.5 
+2.9 
0 
-0.01 


Bus4 
-2.4 
+0.91 
-1.4 
+0.03 
-0.04 
+2.7 
+2.2 
0 
-0.01 


Bus5 
-2.1 
+0.8 
-1.2 
+0.02 
-0.04 
+2.3 
+2 

0 
-0.01 


Bus6 
-2.5 
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-1.4 
+0.03 
-0.05 
+2.8 
+2.4 
0 
-0.01 


Bus8 
-0.6 
+0.33 
-0.2 
+0.12 
+0.11 
+0.3 
+0.3 
+0.01 
0 


3.3. FACTS controllers for system stability enhancement 
Applying SVC and STATCOM at each load bus, one at a time, (3, 4, 5, 6, 8, 9, 10, 19 and 20) 
during fault condition to enhance the stability. For example, the result is taken when fault at bus 8 and 
applying SVC/STACOM at load buses respectively, then the deviation in active power AP and reactive 
power AQ at the load buses respectively due to the fault at bus 8, shown in Figures (8-11) and Tables (6-9). 

From the above data results in the Figures (8-11) and Tables (6-9); during applied of FACTS 
controllers SVC/STACOM, the deviation in active power AP and reactive power AQ at the load buses by the 
fault, is reducing and it becomes near to zero especially when the SVC/STACOM devices applying in bus 9, 
that also prove that it is the optimum location of the FACTS devices, therefore the system returns to stability 
conditions. Also noted from results that the deviation reaches to its minimum and becomes very low incase of 
using STATCOM as compared to using SVC. In other words, STATCOM is better than SVC in terms of 
maintaining the stability of the system. 
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Figure 8. AP in MW due to fault at bus 8 and SVC at load buses respectively. 
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Figure 9. AQ in MVar due to fault at bus 8 and SVC at load buses, respectively 


Table 7. AQ in MVar during fault at bus 8 and SVC at load buses respectively 


VC at AQ at AQ at AQ at AQ at AQ at AQ at AQ at AQ at AQ at 
bus3 bus4 bus5 bus6 bus8 bus9 bus10 bus19 bus20 
Bus3 -1.9 -0.31 -0.3 +0.1 -0.09 +2 +0.4 +0.01 0 
Bus4 -1.2 -4.43 -0.2 +0.11 -0.09 +1.9 +0.4 +0.01 0 
Bus5 -0.8 +0.31 -4.4 -0.03 +0.08 +0.5 +2.5 +0.01 0 
Bus6 -0.9 +0.37 +2.1 +8.69 -0.06 +0.5 +3.9 +0.01 0 
Bus8 -0.8 +0.21 +0.1 +0.09 +10.63 +0.8 +0.9 +0.01 0 
Bus9 -0.5 -0.19 -0.2 +0.1 -0.08 +0.08 +0.4 +0.01 0 
Bus10 -0.8 +0.31 +2.5 -0.17 -0.05 +0.4 +4.2 +0.01 0 
Bus19 -0.7 +0.35 -0.2 +0.12 +0.12 +0.4 +0.3 +0.08 +0.03 
Bus20 -0.7 +0.35 -0.2 +0.12 +0.12 +0.4 +0.3 +0.04 +0.04 
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Figure 10. AP in MW due to fault at bus 8 and STATCOM at load buses, respectively. 


Table 8. AP in MW during fault at bus 8 and STATCOM at load buses, respectivel 





STATCOM AP at AP at AP at AP at AP at AP at AP at AP at AP at 
at bus3 bus4 bus5 bus6 bus8 bus9 bus10 bus19 bus20 
Bus3 +0.8 -0.37 +0.46 +0.19 -0.8 +1.4 -1 +0.1 0 
Bus4 +0.7 -0.16 +0.51 +0.2 +0.8 -1.6 -1.1 +0.1 0 
Bus5 +0.6 -0.41 +0.46 -0.07 +0.7 -1.2 -1.5 +0.1 0 
Bus6 +0.6 -0.41 +0.46 -0.3 +0.6 -1.2 -1.8 +0.1 0 
Bus8 0.6 -0.36 +0.38 +0.19 0 -1.4 -1.3 +0.1 0 
Bus9 +0.07 -0.03 +0.2 +0.1 -0.09 -0.01 -0.01 +0.1 0 
Bus10 +0.6 -0.41 +0.52 -0.08 +0.6 -1.2 -1.8 +0.1 0 
Bus19 +0.6 -0.4 +0.46 +0.21 +0.9 -1.2 -1 +0.3 +0.1 
Bus20 +0.6 -0.41 +0.49 +0.21 +0.9 -1.2 -1.1 +0.2 +0.2 
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Figure 11. AQ in MVar due to fault at bus 8 and STATCOM at load buses respectively. 


Table 9. AQ in MVar during fault at bus 8 and STATCOM at load buses respectively 


AQ at AQ at AQ at AQ at AQ at AQ at AQ at AQ at AQ at 
STATCOM bus3 bus4 bus5 bus6 bus8 bus9 bus10 bus19 bus20 
at 

Bus3 -1.2 -0.06 -0.2 +0.11 -0.1 +1.1 +0.3 +0.01 0 
Bus4 -0.9 -2.18 -0.2 +0.11 -0.1 +1.1 +0.3 +0.01 0 
Bus5 -0.7 +0.3 -2.4 -0.03 +0.1 +0.4 +1.4 +0.01 0 
Bus6 -0.7 +0.31 +1.2 +5.05 -0.08 +0.3 +2.2 +0.01 0 
Bus8 -0.7 +0.25 0 +0.1 +5.67 +0.5 +0.6 +0.01 0 
Bus9 -0.01 -0.12 -0.09 +0.091 -0.07 +0.07 +0.03 +0.01 0 
Bus10 -0.6 +0.29 +1.5 -0.05 +0.08 +0.3 +2.5 +0.01 0 

Bus19 -0.6 +0.31 -0.1 +0.12 +0.12 +0.3 +0.3 +0.03 +0.01 

Bus20 -0.6 +0.33 -0.2 +0.12 +0.12 +0.3 +0.3 +0.02 +0.02 


4. CONCLUSION 

According to above the results, the followings are concluded, Intelligent Fuzzy set technique has an 
important advantage compared with the traditional method indicated as follows: It is an intelligent method 
based on fuzzy rules (if-then statements). It can identify the weakest bus at steady-state conditions and 
transient conditions during high load and fault, respectively, but the traditional method is utilized during fault 
conditions only. From the obtained results, it is clear that Fuzzy technique yields the same results in both 
steady-state and fault conditions which enable the use of the same technique with large systems and applying 
the method in steady-state conditions to get the most vulnerable bus in both steady-state and fault. 
Accordingly, Fuzzy technique is easier, faster and more accurate with increasing the number of 
memberships. Reducing the deviation value for the load buses parameters (P and Q), may reach zero, when 
applying SVC and STATCOM devices respectively at the optimum location (Bus 9 in this paper) during fault 
which means that, the system returned to stability conditions. The degree of reduction in deviation values 
during applying STATCOM, s more than when SVC, that means it is better than SVC for in terms of 
maintaining the stability of the system. When the fault occurs at Bus 20 has more than affected on the system 
because it is the strong bus. 
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